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In condensed matter physics and related areas, topological defects play im- 
portant roles in phase transitions and critical phenomena. Homotopy theory 
facilitates the classification of such topological defects. After a pedagogic in- 
troduction to the mathematical methods involved in topology and homotopy 
theory, the role of the latter in a number of mainly low-dimensional statistical- 
mechanical systems is outlined. Some recent activities in this area are reviewed 
and some possible future directions are discussed. 
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1. Introduction 

Topology is the appropriate mathematical framework for the study of spaces 
which can (and cannot) be continuously deformed into each other. Continuous de- 
formations include twisting and stretching but not tearing or puncturing. Thus a 
cube and a sphere are topologically equivalent entities. Similarly a square is equiva- 
lent to a circle in topological terms. A square is, of course, quite different to a circle in 
that there is a lack of differentiability at its vertices. The appropriate mathematical 
framework to deal with this aspect is differential geometry. Since differentiability 
necessitates continuity, differential geometry is, in a sense, more restrictive than 
topology. However, while the former may yield more concrete results, topological 
arguments generally only lead to existence or classification statements. 

This paper contains a short pedagogic review of certain topological concepts in 
statistical physics, with a focus on homotopy and its consequences. The first part 
(sections EHU) summarizes the basic conceptual and calculational tools used in the 
determination of the homotopy groups for simple topological spaces. The second 
part of the paper (sections EHH1) contains a review of recent progress in statistical 
mechanical models where topological concepts play a crucial role. Conclusions are 
outlined in section fTTJl 



2. Basic notions of topology 

We begin by introducing some basic topological notions, with the primary objec- 
tive of being able to study continuity in mind. We refer the reader to the literature 
(e.g., [ 1,) for basic proofs, which are all rather straightforward. 

Definition 1 If X is a set andT = {X{\ is a collection of finitely or infinitely many 
subsets of X, then we say X is a topological space with a topology T (i) if G T, 
X G T , (ii) if {Xj} is a finite or infinite subset ofT, then the union UjXj G T and 
(Hi) if {Xj} is a finite (not infinite) subset of T , then the intersection HjXj G T. 
The sets X, are called open sets. 

With this definition, the topology consisting only of and X is the one with the 
least number of open sets. For the largest topology, all possible subsets of X and, 
indeed, all points in X are open sets. The latter is called the discrete topology. 
In Euclidean or Cartesian space MJ 1 , one more commonly employs the usual or 
ordinary or Euclidean topology, in which the open sets are restricted to n-balls 
or open intervals. Definition ^ is crucial when analyzing continuity, which is the 
basic purpose of topology. The reason why infinite intersections are not allowed in 
the definition is that such a construct may give rise to open sets under the usual 
topology which consist of a single point. This is not what we usually understand by 
an open interval and would render useless the following definition of continuity. 
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Definition 2 If X and Y are topological spaces and if f : X 
continuous if the inverse of an open set in Y is open in X . 



Y , then f is 



We examine this definition using an example of the usual topology in M, taking 
X = Y = R, and the discontinuous function 



/(*) 



x + 1 if x > , 
x if x < . 




Figure 1. Illustration of the topological 
definition of continuity. 

We now move on to a number of other basic definitions of topology. 



The function is depicted in figure [T] 
It is also demonstrated in the figure that 
while the inverse of some open sets in 
Y are open sets in X, this is not the 
case if the open set in Y is inclusive of 
part of the discontinuous zone. So defini- 
tion 121 captures the notion of discontinu- 
ity. Moreover, a moments consideration 
renders it clear that usage of / rather 
than f~ l in definition |21 would be use- 
less, as / always takes open sets (in X) 
into open sets (in Y). 



Definition 3 A neighborhood of a point x G X is a subset N which contains an 
open set Xi containing e Xi c iv c x. 

Note that N doesn't have to be open (for example, if X = M, with the usual topology 
the closed interval [5,7] is a non-open neighborhood of the point x = 6). But all 
open sets Xi containing x are neighborhoods of x. 

Definition 4 A subset U is closed if its complement X — U is open. The complement 
of an open set is called closed. 

By definition, X is open. Since its complement is 0, the latter is closed. Similarly, 
since by definition is open, its complement X is closed. So and X are examples 
of sets which are both open and closed. 

Definition 5 If U is a set, its closure, written U, is the smallest closed set in which 
U is contained. Because arbitrary intersection of closed sets results in a closed set, 
one may write U = C\ a Vi, where Vi are all closed sets containing U. 

So, for example, the closure of the open set (5, 7) is the closed interval [5, 7]. 

Definition 6 The interior U° of a set U is the union of all open subsets of U . 

Thus the interior of the closed interval [5, 7] is the open set (5, 7) and the interior of 
a closed disk is an open one, for example. 
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Definition 7 The boundary b(U) of a set U is the complement of its interior within 
its closure: b{U) = U — U° . 

Note that U n 6(17) = if U is open. 

Definition 8 A subset U of X is dense in X if its closure U = X . 

For example, with the usual topology on R, there are no open subsets of the set of 
rationals Q. So the interior Q°, which is the union of all open subsets of Q, is 0. In 
this case 6(Q) = Q — Q° = Q. This means that the boundary is the closure. However, 
since b(Q) = R, we have that Q = R. Therefore the set of rational numbers is dense 
in the reals. 

Definition 9 Given a set U and a (finite or infinite) family of sets V = {Vi}, if U 
is contained in UjV; we say that V is a cover of U . If all sets Vi are open then V is 
an open cover. 

For example, the set of open intervals (— n, n), where n G N, is an open cover of R 
under the usual topology. 

Definition 10 A set U is called compact if every open cover has a finite subcover 
(say {Vi, V 2 ,..., V N } ), such that U c V x U V 2 U . . . V N . 

In fact, for a set to be compact, it has to be closed and bounded (and vice versa). 
Essentially this means it has finite volume. 

Definition 11 A set U is called connected if it cannot be written as U = U±{J U 2 
in which Ux H U 2 = 0. It is simply connected or 1- connected if any loop in it can be 
continuously contracted to a point. A domain is a connected open set. 

This definition is illustrated in figure El 
For a set to be simply connected 
it must consist of one component and 
have no holes. Higher-dimensional holes 
are, however, allowed, such as the three- 
dimensional cavity enclosed by the 2- 
sphere, which is a simply connected 
space. 

If the set itself can be continuously 
contracted to a point it is said to be con- 
tractable. All contractable sets are sim- 
ply connected but the converse is not 
true. For example, the sphere is simply 
connected but not contractable - for, if 

it is reduced to a single point, it is no longer a sphere. Contractibility is stronger than 
simply connected and means the space has no holes or cavities of any dimension. 

Our central theme is the study of continuous deformations of spaces from one 
to another. In fact, much information will be gleaned from circumstances in which 




Figure 2. Space X is simply connected 
while Y is connected but not simply (it 
is multiply or non-simply connected) and 
Z is not conected. 
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this is not possible. Indeed, if such continuous deformations are prohibited, then 
there must be some obstacle in the way. These are called topological invariants. 
We attempt such continuous deformations through homeomorphisms ('homeo' and 
'morphism' coming from the Greek words for 'similar' and 'shape', respectively). 

Definition 12 If X and Y are two topological spaces then f : X —> Y is a homeo- 
morphism if it is continuous and if there exists a continuous map g : Y — > X such 
that f o g — ly (the identity map in Y). The map g is also a homeomorphism and 
g o f = l x (the identity in X). I.e., f — g -1 and vice versa. 

With this concept we can categorise 
all topological spaces into equivalence 
classes. Spaces belong to the same class 
if they are homeomorphic to each other. 
To characterize homeomorphic equiva- 
lence classes we need topological invari- 
ants (which are not broken under home- 
omorphism). These include the dimen- 
sion of the space, properties such as 
Figure 3. Representation of a homotopy compactness and connectedness and the 
F(x,t) between two maps f± and fi- powerful concept of homotopy, to which 

we now turn ('homo' and 'topos' come 
from the Greek words for 'same' and 'place'). Homotopy is to continuous maps 
what homeomorphism is to topological spaces - homotopy continuously distorts 
maps while homeomorphism continuously distorts spaces. The concept is illustrated 
schematically in figure 01 in the context of the following definition. 

Definition 13 Let X and Y be topological spaces and f\ and f'2 be continuous maps 
from X to Y . Then f\ is homotopic to f'2 and vice versa if fx can be continuously 
deformed into f2 in the sense that there exists a continuous function F : X x [0, 1] — > 
Y , such that 

F{x,G) = f x {x) , F(x,l) = f 2 (x) . 

Homotopy is an equivalence relation and categorizes all continuous maps from X to 
Y into homotopy equivalence classes which are unchanged under homeomorphism. 

The various so-called homotopy groups provide deep insights into topological 
phenomena in physics, and we now introduce these. 

3. The fundamental or first homotopy group 

If spaces can be continuously deformed into each other without breaking or 
tearing then they belong to the same homeomorphic equivalence class. Clearly, a 
simply connected compact space X (one with no holes) is not in the same homeo- 
morphic equivalence class as a non-simply connected space Y (which has a hole). 
To quantify this rather intuitive statement we consider loops and classes of loops 
within each space. 
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Referring to figure EJ it is clear that 
while loops a and (3 in space X and a' in 
space Y can be shrunk to a point, (3' in Y 
cannot (because the hole encompassed 
by (3' forms an obstruction). We say a 
and 13 are homotopic while a' and j3' are 
not, and write 

a~/3, a' 4 (3' . 







X 



a' 

Y 



Figure 4. All loops in space X are ho- 
motopic. This is not the case in space Y, 
which has a hole. 



In fact every loop in space X can be continuously shrunk to a single point (the 
constant or identity loop) and there is one homotopy class of loops there. In space 
Y, loops can enclose the hole any number neZ times and one may regard clockwise 
orientation as positive n and anticlockwise orientation as negative n. There is an 
infinite number of homotopy classes of loops in Y, each associated with this winding 
number n. 

Loops can be combined (multiplied or added). For example a * (3 is taken to 
mean traversal firstly of loop a and subsequently of loop (3. Inverse loops can also 
be defined, so that a -1 has the same location but reverse orientation to a. The 
product a * a' 1 is then homotopic to (not equal to) the identity loop. 

In fact, since a loop a is homotopic to an infinite number of other loops, it is 
more convenient to consider just one loop, representative of that homotopy class, or 
better, to consider the homotopy class itself - which we label [a] . Products of classes 
are defined as classes of products: [a] * [(3} = [a * /3], and, in particular, [a] * [a' 1 } is 
equal to (not merely homotopic to) the homotopy class of the identity loop. The set 
of homotopy classes defined in this manner has the properties of closure, asociativity, 
the existence of an inverse and an identity. Therefore it forms a group. 

Definition 14 The group of homotopy classes of loops in a topological space X 
based at a point x is denoted by tti(X,x ) and is called the fundamental group 
or first homotopy group. If [a],\J3] G 7Ti(X, Xo), then their product is defined as 
[a] * [f3] = [a * (3] . The identity element is the class of all loops homotopic to the 
degenerate loop comprised soley of the point Xq. 

At this point, we remark that, in the general definition the fundamental group 
is seen to depend on the base point xq. This rather cumbersome burden disappears if 
one restricts one's considerations to pathwise-connected topological spaces. A space 
X is pathwise-connected (also called path-connected or O-connected) if every pair of 
points Xq, Xi G X are connected by a path 7 (i.e., 7 : [0, 1] — > X such that 7(0) = x 
and 7(1) = Xi). This is actually a stronger concept than that of connectedness in 
definition II 11 i.e., connected spaces exist which are not pathwise-connected. For the 
topological spaces encountered herein the two concepts coincide. 

Before stating the theorem establishing the redundancy of a base point in con- 
siderations of the fundamental group structure of most useful topological spaces, we 
recall some basic concepts from group theory. 
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Definition 15 An Abelian group is one in which the elements commute, i.e., if G 
is an Abelian group, for a,b G G, ab = ba. 

For example, Z n is the cyclic group generated by a single element which satisfies 
g n = where is the identity element. All cyclic groups are Abelian. 

Definition 16 A group homomorphism f between two groups G and H is a map 

which preserves the group structure, i.e., for a,b G G, f(ab) = f(a)f(b). 

This definition is sufficient to ensure that the identity is mapped to the identity 
and that the inverse map is preserved. A counter-example in the group of reals 
under multiplication is the sine function, which is not a homomorphism because 
sin (ab) ^ sin a sin b. 

Definition 17 A bijective map f is both infective (also called one-to-one and for 
which a b f(a) ^ f(b)) and surjective (also called onto and for which V6 G H , 
3 a G G s.t. f(a) = b). 

Definition 18 A bijective homomorphism between groups G and H is called an 
isomorphism. We write G = H . 

So a homeomorphism is a continuous isomorphism. (The word 'isomorphism' is 
derived from the Greek, meaning equal shape.) 

Theorem 1 If x ,Xi belong to the pathwise- connected topological space X, then 
■K X (X,X ) = ^i(X,x x ). 

This theorem establishes the independence of tt 1 (X,xq) from xq for reasonable 
topological spaces. We now seek to address the question to what extent does the 
fundamental group depend on X itself. It turns out that it is also independent of X 
up to homeomorphism. I.e., if X is homeomorphic to Y, then 7Ti(X, xq) is essentially 
the same as (isomorphic to) iri(Y,yo). In fact one can do better than that, but 
we need yet another new concept, namely homotopy type, which is broader than 
homeomorphism. It is given by relaxing the equality in definition to homotopy. 

Definition 19 Topological spaces X 
and Y are said to be homotopy equiv- 
alent or of the same homotopy type if 
there exist continuous maps f : X —>Y 
and g : Y — > X such that 

fog~l Y , g o f ~ l x . 

So two homeomorphic spaces are 
necessarily of the same homotopy type, 
but the converse is not the case. Refer 
to figure El where the space Y is formed 




Figure 5. Spaces X and Y (which con- 
sists of a tail T glued on to X at the 
point P) are of the same homotopy type 
but are not homeomorphic. 
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from space X by connecting a 'tail', as illustrated, and let f(x) = x, and g(y) = y 
if y G X, while g(y) — P if y is in the tail. Clearly X and Y are not homeomorphic 
(since / o g ^ ly in the tail), but they are of the same homotopy type. Thus sets of 
the same homotopy type have the same essential structure - they are the same up 
to stretching, twisting or compression but not under cutting. 

Theorem 2 If the pathwise- connected topological spaces X and Y are of the same 
homotopy type, then 

ni(X,x ) = TTi(Y,y ) . 



In particular, if X and Y are homeomorphic, then their fundamental groups are 
isomorphic: iri(X, xq) = 7Ti(Y, yo). Note that the converse does not necessarily hold. 
Nonetheless, this establishes the central result that the fundamental group is a 
topological invariant of a space. 

The fundamental group plays a critical role in the classification of topological 
defects in statistical physics and beyond. Fortunately it is possible to determine this 
group in most cases. To do that, we need the following definition. 

Definition 20 Let R C X. If there exists a continuous map r : X x [0, 1] — > X , 
such that 



r(x, 0) = x V x G X 



r(x, t) = x V x G R 



r(x, 1) G R 




then R is called a deformation retract of X and the map r is called a retract. 

The deformation retract is a sub- 
space of an original space, formed by 
continuous shrinking. In figure the 
space X is a deformation retract of 
space Y. Similarly an annulus, for ex- 
ample can be retracted into a circle. A 
counter example illustrated in figure El 
where the circle C is not a deformation 
retract of the full space X because the 
hole on the right is an obstruction to such a retraction process. The usefullness of 
this lies in the following theorem: 

Theorem 3 If R is a deformation retract of a pathwise connected topological space 
X, then 7Ti(X, xq) — 7Ti(_R, xq). 



Figure 6. The circle C (dashed) is not a 
deformation retract of the full space X. 



Deformation retracts can be used to construct a representative of a space, called a 
polyhedron, which contains all the homotopic features of that space. One can then 
use an algorithm to calculate its fundamental group. This algorithm is illustrated in 
the following examples. 
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Figure 7. Triangulation of an annulus. 



Example 1 (the annulus) To deter- 
mine the fundamental group of an an- 
nulus, one firstly deforms it to a circle 
S 1 and then triangulates that space us- 
ing simplices (see figure EJ). A 0-simplex 
is a point, a 1-simplex is a line segment, 
a 2-simplex is a triangle (including its 
interior), a 3-simplex is a tetrahedron (including its interior), and so on. Here we 
have three O-simplices ({1}, {2} and {3}) and three 1-simplices ({1,2}, {2,3} and 
{3, 1}). The 2-simplex {1, 2, 3} is not involved as the cavity within the triangle is not 
part of the space (the annulus, by definition, has a hole). Each 1-simplex corresponds 
to a group element, which we label g 12 , #23 and g 3 %. Next, construct a 'scaffold' or 
contractable subpolyhedron which contains all the vertices of the triangulation or 
polyhedron. In our case {1, 2} U {2, 3} will suffice. 

This contractable subpolyhedron is 
illustrated for the present example in 
figure |H1 Each of the 1-simplices in the 
spanning subpolyhedron gives the iden- 
1 tity group element which we denote by 

0, so gi2 = #23 = 0. We are left with 
Figure 8. A subpolyhedron (dashed) one non-trivial element, namely g 3 i = g. 
spanning the triangulation of S 1 . The group generated by such an element 

is Z, the integers under addition. There- 
fore the fundamental group of the annulus is vri(S' 1 ) = Z. 



Example 2 (the disc) The disc D or D 2 is again triangulated by a triangle, but 
this time its interior is included. This is the 2-simplex {1,2,3}. Such a 2-simplex 
gives a relation gi2g22,9^i — 0, which, with g 12 = #23 = gives g 3 i = 0. So all 
elements of the group are the identity and H\{D) = {0}. 





Example 3 (the Mobius strip) The 

Mobius strip, a suitable triangulation 
and a contractable subpolyhedron are 
illustrated in figure M Each of the 
five 1-simplices contained in the sub- 
polyhedron give the identity group ele- 
ment. There are seven remaining group 
elements corresponding to the seven 
1-simplices outside the subpolyhedron 
(taking care to note the left and right 
edges are the same). But there are six 2-simplices giving six relations between these 
seven elements. Only one non-trivial element remains and m (Mobius) = Z. 

For any of these simple spaces, we can guess the fundamental group by imagining 
whether or not loops straddling the structure can be deformed to a point. 



Figure 9. The Mobius strip (left), its tri- 
angulation (right) and a spanning sub- 
polyhedron (dashed). 
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Example 4 (the projective plane) 

The real projective plane RP 2 is con- 
structed by identifying each pair of dia- 
metrically opposed points on the bound- 
ary of a disc (figure ITn|l . A path connect- 
ing a point to its antipode is therefore Figure 10. Formation of the projective 
closed and constitutes a loop which can- piane RP 2 by identification of antipodal 
not be shrunk to a point. A trivial loop points on a disc, 
is only constructed by returning to the 

starting point (not antipode). So, besides the identity, there is a nontrivial group 
element, which is its own inverse and vr^RP 2 ) = Z 2 = {0,g}. 

To determine the fundamental group of higher-dimensional spaces, van Kampen's 
theorem is useful: 

Theorem 4 If X andY are topological spaces, the fundamental group of their prod- 
uct is the direct product of their fundamental groups. 

miX x Y) ^7ri(X)©7ri(y) . 

For example, the three-dimensional torus is given by T 2 = S 1 x S 1 so 71"! (T 2 ) = 
MS 1 ) ©^(S 1 ) = ZffiZ = Z 2 . 




4. The higher homotopy groups 

Only 1- and 2-simplices are used in the algorithm to calculate the fundamental 
(or first homotopy) group. This means that 7Tl can only detect two-dimensional 
holes present in a space (or combinations thereof, such as the two two-dimensional 
holes in the torus). For example, by considering shrinkage of a loop on its surface, 
it is easy to convince oneself that the fundamental group of a sphere in 3-space is 
trivial: 7r 1 (S' 2 ) = {0}. Such one-dimensional loops are incapable of detecting the 
three-dimensional interior of the sphere. To deal with this circumstance we need 
higher homotopy groups. 

While the reader is referred to the literature for the mathematial details [Q], the 
basic idea is straightforward and now outlined. The 1-loops we encountered hitherto 
can be thought of as elastic bands tied at a basepoint xq. The two-dimensional 
equivalent of such an object can be thought of as a balloon, with the neck anchored 
at xq. Such a 2-loop can encompass a spherical hole in much the same was as a 1-loop 
can enclose a circular one. Moreover, a balloon can be wrapped or unwrapped an 
integer number of times around a sphere and a group structure can be built up in 
an equivalent manner to the (homotopy equivalent classes of) 1-loops. In still higher 
dimensions, n-loops can be defined. 

Definition 21 The homotopy class of n-loops in a topological space X with base- 
point xq forms the n-dimensional homotopy group 7r n (A, Xq). 
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Table 1. The homotopy groups of the n-sphere. The higher groups are Abelian 
but fall into no easily distinguishable pattern. 



m 


n 

1 


2 


3 


4 


5 


6 


7 


8 ... 





{0} 


{0} 


{0} 


{0} 


{0} 


{0} 


{0} 


{0} ... 


1 


z 


{0} 


{0} 


{0} 


{0} 


{0} 


{0} 


{0} ... 


2 


{0} 


z 


{0} 


{0} 


{0} 


{0} 


{0} 


{0} ... 


3 


{0} 


z 


z 


{0} 


{0} 


{0} 


{0} 


{0} ... 


4 


{0} 


z 2 


z 2 


z 


{0} 


{0} 


{0} 


{0} ... 


5 


{0} 


z 2 


z 2 


z 2 


z 


{0} 


{0} 


{0} ... 


6 


{0} 


Zl 2 


Zl2 


z 2 


z 2 


z 


{0} 


{0} ... 


7 


{0} 


z 2 


z 2 


z©z 12 


z 2 


z 2 


z 


{0} ... 


8 


{0} 


z 2 


z 2 


z 2 ©z 2 


Z 2 4 


z 2 


z 2 


z ... 


9 


{0} 


z 3 


z 3 


z 2 ©z 2 


z 2 


z 24 


z 2 


z 2 ... 



It turns out that, while some but not all fundamental groups are Abelian, all higher 
homotopy groups are Abelian. While there is no algorithm to calculate the higher 
homotopy groups comparable to the case of the fundamental group, nor is there 
a higher homotopy analogue to van Kampen's theorem, there are tools such as 
Freudenthal's theorem: 

Theorem 5 The m- dimensional homotopy group of the n-sphere ir m (S n ) depends 
only on m — n for m < 2(n — 1) . 

An immediate consequences of this is ir m (S n ) = n m+ i(S n+1 ) for m < 2(n — 1). 
Further, since 7Ti(5 n ) = {0} for n > 2, 7c m (S n ) = {0} for n > m+ 1. If n — m, one 
has (for n > 2) ir n (S n ) S vr 2 (S 2 ) S Z. 

Finally, although tt is not a group, it is often used to denote the number of 
connected domains or components in a space. 

Classification of the homotopy groups of topological spaces is an active field 
of research and has generated many surprising results. Table [0 illustrates this by 
listing Tr m (S n ) for small m and n. While there is no overall identifiable pattern, as 
the dimension increases some regularity does occur. In particular, 7r n+ i(5' n ) = Z 2 
for n > 3, TT n+2 (S n ) = Z 2 for n > 2, n n+?J (S n ) = Z 24 for n > 5 and so on. 

Homotopy groups for other topological spaces have also been determined and 
some of those more commonly used in physics are listed in table 121 

This completes the pedagogic introduction to homotopy. Equipped with the 
above topological notions, one may now examine the role played by homotopy in 
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Table 2. The homotopy groups of various useful topological spaces. 



7Ti(R*) = 7n(D) = 7T m (5"> m ) VT^S 1 ) {0} 

tt^T 2 ) = Z©Z = Z 2 

71" i (T n ) = Z n 

ttx (RP n ) = Z 2 (n > 2) 
vr m (RP n ) S vr m (5 n ) for m > 2 

7T rt (RP n ) = Z 



statistical physics. 



5. Phase transitions and topological defects 

Phase transitions are amongst the most ubiquitous and remarkable phenomena in 
nature and involve abrupt or gradual changes in quantifiable macroscopic properties 
brought on by varying a system's parameters such as temperature, pressure or a 
coupling. Equilibrium statistical physics which describes such phenomena is based 
on the premise that the probability that a system is in a state S with energy E at 
a temperature T is 

= (1) 

where (3 = 1/ksT, ks is a universal constant, and Zl(/3) is the partition function 
which serves as a normalising factor. It is given by 

Z L ((3)=J2exp(-f3E(S)). (2) 

s 

Here, and in the following, the linear extent of the system is indicated by the sub- 
script L. Another fundamental quantity is the free energy, given by 

f L ({3) = ±- d ln Z L (f3), (3) 

where d is the dimensionality of the system so that its volume is L d . In the modern 
classification scheme, phase transitions are categorised as first, second or higher 
order if the lowest derivative of the infinite- volume free energy foo(P) that displays 
non-analytic behaviour is the first, second or higher one (see, e.g., [ 2J). 

For an 0(n)-symmetric model, the sites i of a <i-dimensional lattice are occupied 
by n-dimensional unit-length vectors s*j G 5* n_1 which may be considered to reside in 
an internal spin space. In general then, n may be independent of the dimensionality 
d of the physical space or medium. For the 0(n) model, the energy of a given 
configuration is given as 

E = -J2 SiSj , (4) 

(ij> 
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where the summation runs over nearest neighbouring sites or links of the lattice. 

The n — 1 version of this construct is the Ising model. If n = 2, the spins live in 
a plane and the model is referred to as the XY model. The 0(3) version is called the 
Heisenberg model and the limit n — > oo corresponds to the spherical model. In each 
of these models, the system (and the partition function in particular) is invariant 
under rotations in spin space. However, the order parameter or spin expectation 
value (sj) may, in principle, not respect this symmetry. This very common circum- 
stance is referred to as spontaneous symmetry breaking. Let T c or f3 c denote the 
phase transition point. The conventional scaling scenario at such a transition is of 
the power-law type in the reduced temperature t = T/T c —1. For example, this is the 
type of scaling which characterises the standard temperature-driven paramagnetic- 
ferromagnetic phase transition in the XY model in d > 2 dimensions. 

However, the existence of a phase with long-range order is precluded in two- 
dimensional models with continuous symmetry group and continuous interaction 
Hamiltonian, such as (JH), due to the Mermin- Wagner theorem [Ej. Nonetheless, a 
temperature-driven transition can still exist - to a phase with topological order. 
The two-dimensional XY model has a famous infinite-order phase transition of this 
type which breaks no system symmetries and is called the Berezinskii-Kosterlitz- 
Thouless (BKT) transition [IHE]. The transition is mediated by topological defects 
called vortices. Such vortices can be identified by tracking the well-defined values 
of the spins Sj along a given contour (see figure [TTJ) . Spins can rotate along such a 
contour through 2mr for n e Z. If n ^ 0, one speaks of a vortex (actually one may 
refer to a negative-^ vortex as an antivortex). In the absence of the lattice, shrinking 
the contour to a point would lead to a singularity in the presence of such a vortex. 

From homotopy theory, it is clear that the winding number in the above two- 
dimensional example is related to the fundamental group of the order-parameter 
space 7r 1 (S' 1 ). In principle, in a d = 2 dimensional medium, one could also have a line 

defect if tt was non-trivial (since ttq 
denotes the number of connected do- 
mains in the medium, such a defect 
is called a domain wall). In a d — 3 
medium, defects of dimension 0, 1 or 
2 may exist. These are point defects 
(monopoles), line defects (vortices) or 
wall defects (domain walls), and are de- 
tected by 7r 2 , 7Ti and 7T respectively. In 
general, then, an m-dimensional defect 
in a d dimensional medium is classified 
by ^d-m-i(S), where S is the order- 
parameter space. They are called do- 
main walls (classified by 7r ), vortices 
(tTi), monopoles (n 2 ), instantons (7r 3 ), 
etc. The classification scheme is summa- 
rized in table 3. 
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Figure 11. Examples of an antivortex (or 
vortex of charge —1) and a vortex (of 
charge 2) in the XY model on a regu- 
lar lattice. As contours around the two 
bold points are traversed (counterclock- 
wise), spins rotate through — 2n and Air 
respectively. 
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Table 3. Classification of topological defects. Here d is the dimensionality of the 
medium and m is that of the defect. 
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6. The BKT transition in the XY and step models 

The XY model in two dimensions is the generic one for the study of the topo- 
logically-mediated BKT transition. It is also used to study systems such as films 
of superfluid helium, Josephson-j unctions, superconducting materials, fluctuating 
surfaces as well as certain magnetic, gaseous and liquid- crystal systems. Besides 
this model, transitions of the BKT type also exist in the ice-type F model [E], 
antiferromagnetic models [Ej, certain models with long-range interactions [|H], lattice 
gauge theories [ Ej and in string theory [ EJj amongst others. Thus a thorough 
and quantitative understanding of the paradigmatic two-dimensional XY model is 
beneficial to a number of areas within theoretical physics. The BKT transition in 
this model has come under intense scrutiny, and only in recent years has the picture 
begun to become clear. 

At low temperature, vortices and antivortices appear in pairs and there is quasi- 
long-range order (meaning orientation of spins over small length scales). The corre- 
lation function characterises this ordering and, with Xj denoting the position of the 
i th site, its leading critical behaviour from renormalization group (RG) theory is 

Goo(x 4 - x 3 ) = (s iSs ) ~ \ Xi - X] \- {d ' 2+ ^ (5) 

(with d = 2 in this instance). The correlation length in an infinitely large system is 
denoted by £oo(t) and measures to what extent spins at different sites are correlated. 
It diverges in the massless low-temperature phase. This divergence persists, with the 
system remaining critical with varying rj(/3), up to f3 = (3 C at which r](/3 c ) = rj c = 1/4. 

According to the BKT picture, the number of vortex pairs increases as the tem- 
perature is raised and they proliferate in the high-temperature phase. At the transi- 
tion point, it becomes no longer sensible to describe vortices as belonging to paired 
sets and they are said to unbind. In the high-temperature phase, the dissociation 
of vortices destroys the order of the system. Above this point, correlations decay 
exponentially fast with leading behaviour 

Goo(xj - Xj) ~ exp [-\xi - Xj\/£oo(t)] . (6) 
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The leading approach to the transition from the high-temperature phase is charac- 
terised by essential singularities, 





~ exp ( 




t-2 
Soo 1 


Xoo(f) 


Soo 



for the correlation length, the specific heat and the susceptibility (which respectively 
measure the response of the system to variations in the temperature and application 
of an external magnetic field). Here b is a non-universal constant. 

The 'BKT scenario' is thus taken to mean a phase transition which (i) exhibits 
essential scaling behaviour and (ii) is mediated by a vortex-binding mechanism. This 
picture is based on perturbation theory [HJE] and a variety of techniques have been 
employed for over thirty years in attempts to verify it. 



6.1. BKT scaling behaviour 

These attempts at verification have included non-perturbative simulational and 
high-temperature series analyses of the XY model. In particular, verification of the 
analytical BKT RG prediction that v — 1/2 and r] c = 1/4 and accurate deter- 
mination of the value of f3 c proved elusive. Typically the critical temperature was 
determined by firstly fixing v = 1/2. However, measurements of r\ c then yielded a 
value incompatible with the BKT prediction. The elusiveness of such unambiguous 
corroborative evidence led to the essential nature of the transition in the d = 2 XY 
model being questioned [ 111 j . An extensive overview of the status of the model up 
to 1997 is contained in [I12j. 

A self-consistency analysis, based on Lee- Yang zeros, was used in [ [TJl [T3] to 
show that the thermal scaling forms (jHJ) and (and similar formulae for other 
quantities) are mutually incongruent as they stand, and have to be modified to 
include multiplicative logarithmic corrections. In fact, 

CUO ~ U(t)- 2 (\nU(t))\ (10) 

Xoo(t) ~ Uitf^^Uit))- 2 ' , (ii) 

and 

(j qq £ 3i * j J I X <i I \ | *^ i *^ j I y i ( ^~ ^ ) 

with q = 6 and where RG indications implicit in are that r = —1/16 = —0.0625. 
Besides this implication, the first indication of the existence of a non-zero correction 
exponent r was presented by Butera and Comi using high-temperature series ex- 
pansions [El- On the other hand, arguments that BKT theory implies the absence 
of multiplicative logarithmic corrections, i.e. r = 0, have also been recently made [ 

UniEEE]- 

The formulae (jl(Jj) and (jlip refer to thermal scaling on infinite lattices. Such systems 
are not amenable to computational techniques as one may only simulate finite- 
size systems. There, finite-size scaling (FSS) theory predicts that the role of the 



15 



R. Kenna 



Table 4. Estimates for the critical point (3 C for the XY model as well as the 
logarithmic-correction exponent r for the BKT universality class from a selection 
of recent papers. 



Authors 


Year 


Method 


Pe 

f t- 


T 


JCostprlitz TViouIpss flnl 


1973 


RG 




— 0625 


Irving k Kenna [H2HH1 


1996 


FSS 


1.113(6) 


-0.02(1) 


Patrascioiu & Seiler [ I21| 


1996 


thermal 




0.077(46) 


C^rnmnostrini Pelissetto 


1996 


t, Tiprmal 

Lf AAV^-i. lllCil 


1 1158(6) 


042(5) 


Rossi & Vicari [HI] 






1.120(4) 


0.05(2) 


Janke [HTj 


1997 


FSS, 




-0.027(1) 






thermal 




0.0560(17) 


Hasenbusch & Pinn [ 23 


1997 


RG 


1.1199(1) 




Jaster & Hahn [ITS] 


1999 


FSS, 




-0.0233(10) 






thermal 




0.056(9)/0.070(5) 


Balog, Niedermaier, Niedermayer, 


2001 


RG, 







Patrascioiu, Seiler, Weisz [[ToUTo] 




thermal 






Tomita & Okabe [EH 


2002 


FSS 


1.1194(8) 


0.038(5) 


Dukovski, Machta & Chayes [05] 


2002 


FSS 


1.120(1) 




Chandrasekharan & Strouthas [I19j 


2003 


FSS 




-0.035(10) 


Hasenbusch [ 26 


2005 


FSS 




-0.056(7) 



correlation length is played by the lattice extent L. For example, the susceptibility 
at criticality (t = 0) scales as [ IT5] 

XL (0)~L 2 -^(\nL)- 2r . (13) 

Verification of the BKT scaling scenario then requires numerical confirmation of 
the thermal formula (jllj) and/or the FSS formula (|13|) . In [ [121 E], the hitherto 
conflicting results for u, (3 C and r] c were finally resolved. However, the analysis re- 
sulted in an estimate of —0.02(1) for r, a value in conflict with the RG prediction 
of r = -1/16 = -0.0625 from [EJ. 

This prompted the shifting of the focus of recent numerical studies of the XY 
model to the determination of the logarithm exponent r. Indeed, in [[TTJ [TH1 [T^] , FSS 
analyses using (fl3|). yielded values compatible with that of [H21IIB] but incompatible 
with [[3] (see table 4 and [ 20 for an overview). Nonetheless, it was clear that the 
resolution of the u-f3 c -r] c controversy is achieved by taking the logarithmic corrections 
into account. The most precise numerical estimate for the critical temperature for 
the XY model is (3 C = 1.1199(1) [I23] 5 a value obtained by mapping the model 
onto the body-centered solid-on-solid model (which is exactly solvable), and thereby 
circumventing the issue of logarithmic corrections. Indeed, the recent analyses which 
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have included logarithmic corrections have resulted in estimates for (3 C compatible 
with this value and with rj — 1/4 for the anomalous dimension. 

The study of the XY model contained in [ E] employs the Villain formulation, 
which has an action different to (@J) and leads to a different value for quantity f3 c 
(which is nonuniversal) . Similarly, in table 4, [ and contain studies of a 
lattice grain boundary model and a lattice gauge theory respectively. Although the 
critical temperatures for these models bear no relationship to their XY counterpart, 
they have the same scaling behaviour as they are in the same universality class. Thus 
the value for the logarithmic exponent r should be the same as in the XY model. 

It was further suggested in [^JE] that since taking the logarithmic corrections 
into account leads to the resolution of the leading-scaling puzzle, in the same spirit 
numerical measurements for r may become compatible with the RG prediction if 
sub-leading scaling corrections are taken into account. In this case, (j!3|) is more 
fully expressed as 

X,(0)~L 2 -MlnLr-|l + O^)} . (14) 

This was tested in [ H21 HZj, but the lattice sizes available (up to L = 256) were 
too small to resolve the issue. (For a similar situation in the q = 4 two-dimensional 
Potts model, see [I2ZI-) The problem was revisited recently in [|2Sj in which a very 
high precision numerical simulation achieved lattices as big as L = 2048. Using the 
alternative ansatz 

X L(0)~L 2 -^(C + \nL)- 2r , (15) 

with C an additional free parameter, this analysis gave r = —0.056(7) provided 
suitably large lattices are used (with L > 256). This value is compatible with the 
analytic prediction r = —0.0625. In [Ej, a very careful study of the F model, which 
is exactly solvable and has a BKT-type phase transition, further elucidates the 
necessity using sufficiently large lattices when numerically analysing models with 
highly subtle logarithmic corrections. 

It is interesting to note that, with the exception of [121], all of the FSS-based 
analyses of table 4 yield a negative value for r, in line with the BKT prediction. It 
is rather surprising that, in table 4, all of the thermal scaling analyses (as well as [ 
24J) yield positive values of r - far from the RG prediction that r = —0.0625. This 
indicates that the FSS approach may generally be more powerful than the thermal 
scaling one, although more extensive computational analyses would be required to 
compare this approach to 

Recently, Berche et al. have used a conformal map to rescale distances on the 
lattice so that it is less sensitive to finite-size boundary effects and then to determine 
T}(0) at any temperature in the critical phase [|2IIJI2Ell2n]- In particular, their numer- 
ics yield accurate agreement with the analytic prediction 1] C = 1/4 at the transition 
point. It is interesting to note that this agreement was achieved very straightfor- 
wardly using moderate computational effort and without recourse to enormous lat- 
tices. Rather surprisingly, the evidence showed no noticable logarithmic corrections 
to the order-parameter density profile at the transition [I2E1I2H] (see also [ED])) but 
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was clearly supportive of the existence of such corrections in the correlation func- 
tion EH]- However, the lattices used were not sufficiently large to unambiguously 
confirm the quantative behaviour of the logarithmic corrections. 

6.2. BKT vortex unbinding 

The second crucial aspect in the BKT scenario is the vortex unbinding mech- 
anism. For the XY model the energy of a single isolated charge-n vortex on an 
L x L lattice is proportional to n 2 InL/a where a is the lattice spacing (which never 
vanishes in a real physical system made up of atoms). The total energy of two 
vortices of charge n and — n centred at x\ and £2 is, however, proportional only 
to n 2 ln|xi —X2\ja. At low temperatures, then, vortices occur mostly in vortex- 
antivortex (dipole) pairs. Mutual cancellation of their individual ordering effects 
means that such a pair can only affect nearby spins and cannot significantly dis- 
order the entire system. Topological long-range order exists in the system at low 
temperature. As the temperature is raised, then, the vortices proliferate and the 
distance between erstwhile partners becomes so large that they are effectively free. 
I.e., the BKT transition is one from a phase of dipoles to a plasma of vortices, which 
render the system disordered. 

It was long believed that adjusting the vortex-binding dynamics of the XY model 
may disable the BKT transition mechanism, leading to a different one or the absence 
of a transition of any type [ EH] . With this in mind, the step model is derived from 
the XY model by replacing the Hamiltonian (jlj) with 

E = - sgn (siSj) . (16) 

While the configuration space for this model is globally and continuously symmetric, 
as in the XY case, its interaction function is discontinuous. The energy associated 
with a single vortex for this system is expected to be independent of the lattice 
extent, leading to the expectation that the disordered vortex-plasma phase may 
exist for all temperatures. If this is the case, there was expected to be no vortex- 
driven phase transition in the model; if there is a phase transition, it was expected 
not to be of the BKT type. Indeed, early studies supported this assertion (see [ IT2] 
for a review up to 1997). 

In [ ^] , strong and clear evidence was presented that a transition exists in 
the model, and, that it belongs to the same universality class as the XY model 
in two dimensions (with even the corrections to scaling, insofar as they could be 
discerned, being the same as those for the XY model). Because of the dissimilar 
vortex energetics of the two models, this came as a surprise. The question then arose 
as to how a transition mediated by vortices (in the XY model) can be insensitive 
to the energetics of such vortices. 

In the Mermin- Wagner theorem was extended to discontinuous interaction 
functions including that of the step model The issue of the phase transition 
there was again addressed in [ [33] , where further evidence for the existence of a 
BKT transition in the step model was presented. In the approach was to focus 
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on numerical analyses of the helicity modulus, which experiences a jump at the 
transition (for a similar approach to the XY model see [El])- The main point of [ 
13*3*] is that the cost in free energy for a single vortex is proportional to \nL/a and 
that this is the feature which stabilizes the low-temperature (dipole) phase. I.e., this 
inhibits proliferation of free vortices in the low-temperature phase. It further explains 
the occurence of a transition in the step model. The support for the existence of 
a BKT transition in the step model given in [ EE] indicates that the BKT vortex 
scenario is more general than was heretofore realized. 

7. The diluted XY model 

A vibrant current area of research is the question of the role and consequences 
of the presence of impurities in various systems, including the XY model. The 
occurrence of physical impurities renders any model more realistic, as such defects 
are often present in actual (and porous) systems. These physical impurities are 
modelled by removing (diluting) the sites or bonds of the lattice. Clearly, if the 
dilution is sufficiently strong the percolation of spin-spin interactions across the 
lattice is curtailed and it is effectively broken into finite disconnected sets. This 
occurs at the percolation threshold. In that circumstance, no phase transition can 
occur for any model, as a true phase transition necessitates a thermodynamic limit. 
More moderate dilution is generally expected to lower the transition temperature 
relative to its value for a pure (undiluted) system. 

The special feature of the XY model is the presence of vortices as the mecha- 
nism mediating the transition. It turns out that vortices are attracted to the physical 
impurities (vacant sites or bonds) and the vortex energy is reduced at such a va- 
cancy. Therefore as the dilution is increased, more vortices can be formed and the 
consequent amount of vortex-induced disorder in the system is increased. This in 
turn may enhance the lowering of the critical temperature to such an extent that it 
vanishes even before the percolation threshold is reached. 

This issue is addressed in [ EI3 EH1 EH EE]- According to the Harris criterion, 
disorder does not change the leading scaling behaviour of a model if the critical ex- 
ponent a associated with the specific heat of the pure model is negative [El] - This is 
the case for the XY model in two dimensions [ES]- Thus the critical temperature in 
the diluted d = 2 XY model can be safely (albeit approximately) identified as being 
the location at which r](/3) = 1/4. The vanishing of the critical temperature then 
gives the critical vacancy density in that model. In [ EE] the critical temperature 
was reported to vanish at site- vacancy density p ~ 0.3. However, for d = 2 regular 
lattices, the percolation threshold (the density of vacancies required to disconnect 
the lattice) is p ~ 0.41. Thus the work of [ EH] suggests that, indeed, the critical 
temperature vanishes before the percolation threshold is reached. In a later study, 
however, Berche et al. suggested that the critical site-density is closer to the perco- 
lation threshold [EZ]- Support for the latter result also recently appeared in [EE]- 
These more recent studies - which concern the site-diluted version of the model - 
indicate that the vortices do not, in fact, strongly enhance the lowering of the critical 
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temperature in the d = 2 XY model. 

This conclusion is further supported by a similar recent study [HU] for the bond- 
diluted XY model, which favours the general scenario depicted in [|37J|3B! over that 

of [Eg. 

The above studies concern the value of the critical temperature in diluted two- 
dimensional XY models. The precise scaling behaviour of the thermodynamic func- 
tions at the phase transition is also of interest. The Harris criterion concerns only 
the leading scaling behaviour (which does not alter in the XY models considered 
here) and does not predict what effect dilution can have on the quantitative nature 
of the exponents of the logarithmic corrections in the negative-a case. This would 
be an interesting avenue for research in the future, and the two-dimensional XY 
model offers an ideal platform upon which to base such pursuits [I41j. 

8. 0(n) and RP n 1 models 

In this section, a selection of recent results in non-Abelian 0(n) models in two 
dimensions and in RP n_1 models in two and three dimensions are discussed, focusing 
on certain aspects that are still unresolved. For work on the 0(n) models in three 
(as well as two) dimensions, the reader is referred to the literature [ I42| . 

8.1. Non-Abelian 0(n) models in two dimensions 

It is widely believed that there are fundamental differences between models with 
Abelian and non-Abelian symmetry groups. The 0(2) symmetry group of the XY 
model is Abelian and all 0(n) groups with n > 2 are non-Abelian. From the Mermin- 
Wagner theorem, any continuous symmetry of the 0(n) type cannot be broken in 
two dimensions [Ej. On this basis, there cannot be a transition to a phase with 
long-range order in any 0(n) model with n > 2 there. (The n = 1 case is discrete, 
and the corresponding Ising model possesses an ordered phase in d = 2 [I43j.) As we 
have seen, in a two-dimensional theory, topological defects of dimension m can exist 
if the (1 — m) th homotopy group, vri__ m , of the order-parameter space (which for 0(n) 
models is the hypersphere S 1 ™ -1 ) is non-trivial. From table ^ the only non-trivial 
homotopy group of the form 7r 1 _ m (S' n ~ 1 ) is vr 1 (5' 1 ) = Z. This is the condition that 
gives rise to point defects (vortices) with integer charge in the n = 2 case (the d = 2 
XY model). The binding of these vortices at low temperature is the mechanism 
giving rise to the BKT phase transition 

For n > 2, therefore, conditions are not supportive of the existence of stable 
topological defects of this type in two dimensions and the majority belief is that 
there is no distinct low-temperature phase and consequently no positive-temperature 
phase transition in these models [HI]. There is a vast literature on the d = 2 0(n) 
models with n > 2 and the reader is referred to [ H2] for a review. Furthermore, 
while perturbation theory predicts that these n > 2 models are asymptotically 
free, there is no rigorous proof to this effect. (Asymptotic freedom means that the 
effective strength of interactions vanishes as the energy is increased.) This notion has 
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been questioned and in [ [TBI EJJ HE] evidence in support of the existence of phase 
transitions of the BKT type in these models has been given, as well as heuristic 
explanations of why such transitions could occur and a rigorous proof that this 
would be incompatible with asymptotic freedom. 

Perturbative and high-temperature series expansions [ JJ] as well as Monte Carlo 
calculations for n > 3 [HB] have been performed, and these do not support the exis- 
tence of such a BKT-like phase transition in these two-dimensional models. Instead 
they are in agreement with perturbation theory and the asymptotic freedom sce- 
nario. Nonetheless, the controversy has not been entirely resolved [ |17j and it is 
plausible that inclusion of logarithmic considerations in these considerations could 
help in the search for a precise unambiguous resolution. 

8.2. Liquid crystals and RP n_1 models 

The liquid-crystal state is a phase of matter which exists distinct from (i.e., not a 
mixture of) the solid and liquid states. Unlike in a normal liquid, which is isotropic, 
in a liquid crystal the properties are directional dependent. This is because the 
molecules of the liquid crystal are elongated in one direction in three-dimensional 
physical space. On the lattice, each molecule may be represented by a directed rigid 
rod. Such a direction is without orientation, so that the system is unchanged by 
flipping the rod. For the n-component model the Hamiltonian is 

E = -j2m)\ (17) 

(id) 

where Si are n-component unit vectors and the quadratic form gives that — s« and 
Si describe the same direction [EE]- The symmetry group for this n- vector model is 
RP n -\ This is the space in n dimensions formed by identifying antipodal points on 
an (n — l)-sphere S 1 " -1 . Alternatively, restricting to one of the hemispheres of S 1 " -1 , 
the real projective space RP n_1 can be considered as equivalent to the (n — l)-disk, 
D n ~ l , with antipodes on the boundary S n ~ 2 identified (see example 4 of section EJ). 

For n = 2, the real projective space RP 1 is topologically equivalent to the circle 
S 1 in 2-space and the trigonometric identity 2 cos 2 9 = 1 + cos (26) leads to the model 
(f*T7j) being equivalent to the XY model. 

In three dimensions, then, the condition for the existence of topological defects 
in the n-component model is the non-triviality of 7T2_ m (RP n ~ 1 ). Since RP n_1 is 
connected and 7To(RP n_1 ) is trivial, there are no two-dimensional (wall) defects in 
a three-dimensional medium for any n. In the liquid crystal case where n = d = 3, 
table |21 gives that 7r 2 (RP 2 ) = Z so that point defects (monopoles, called hedgehogs 
in this instance) may exist. (In addition to hedgehogs, if the system is finite in extent, 
there may exist pointlike topological entities called boojums, which loosely speaking 
are like half- hedgehogs - see [HI] and references therein.) Also from table 121 since 
7Ti(RP 2 ) = Z 2 , line defects (called disclinations) can exist in the n — d — 3 system. 
Indeed, this is reflected in the tendency of such molecules to align themselves into 
linear or threadlike patterns. Such materials, which can cause polarization of light, 
are called nematic liquid crystals after the Greek prefix 'nemato' meaning threadlike. 
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Actually there are other types of liquid crystal states beyond the nematic one and 
the reader is refered to the (vast) literature [150]. 

Similarly, in two dimensions, point defects (vortices) owe their existence to the 
non-triviality of 7ri(RP n_1 ), which is isomorphic to Z2. Even though only point 
defects can exist in the two-dimensional version, one often also generically refers to 
a nematic phase here too. 

Perturbative RG analyses of the RP n_1 model (which neglect the effects of topo- 
logical defects) predict no transitions in two dimensions [|Hj and a second-order one 
in three dimensions for n > 3 [EIJ. However, numerical and experimental research 
indicates that the three-dimensional RP 2 transition is, in fact, weakly first-order [ 
HB1 E21 EH] • The three-dimensional RP 3 model also has a first-order transition [ 151] 
and is closely linked to frustrated spin systems which have been extensively studied 
(for a review of theoretical, numerical and experimental work, see [I55j). The large 
n limit of the d = 3 RP n_1 model also possesses a first-order transition [EH]- It 
seems clear that topological defects are responsible for the first-order phase tran- 
sition in both the three-dimensional RP n_1 models and their frustrated magnetic 
counterparts [155]. 

Since the RG approach fails in the three-dimensional model (where there is, in 
fact, a first-order transition), it is possible that it also fails in two dimensions and 
that this failure is due to the transitions being topologically driven [|57j. Indeed, it 
has been shown that, in two dimensions, while perturbative RG predictions match 
well with Monte Carlo measurements when the model has trivial homotopy, there 
is a clear disagreement between the two approaches when topological defects are 
present [155]. 

The two-dimensional RP n_1 models were considered in [|57], with n = 3 and 
n = 40. In the nematic n = 3 case, evidence was presented for a transition described 
by a diverging correlation length and susceptibility but a cusp (as opposed to a di- 
vergence) in the specific heat was reported. Similar to the two-dimensional XY case, 
both the correlation length and the susceptibility appeared to remain infinite below 
the critical temperature. Despite the scale of the study performed in [|57j, it was 
not possible to distinguish between essential scaling and standard power-law scaling 
as the transition is approached from the high-temperature phase. The importance 
of the defects was, however, clearly demonstrated, in that they carry most of the 
energy and the transition appears to be again mediated by their unbinding with 
increasing temperature. The d — 2, n — 3 case was also analysed in [.59], where 
a second-order phase transition was favoured. In [ |57] , a numerical analysis of the 
n = 40 case was also compared with analytical results for n = 00. The latter has 
a topologically mediated first-order transition for all d > 2 [ 56J. The possibility of 
(lattice-dependent) first-order transitions at large or infinite N was discussed in [ 
ED]- 

The powerful conformal techniques used in [1201123121] were similarly employed 
in [EI]) favouring a nematic/isotropic topologically-mediated transition in the two- 
dimensional RP 2 model and a close similarity with the two-dimensional XY model. 
The effect of the suppression of the topological defects was explored in [E2] where it 



22 



Homotopy in statistical physics 



Table 5. Summary of the status of standard 0(n) and RP n 1 models and some 
recent developments considered in this paper. Here m is the defect dimension. 



Model 


d 


Homotopy group, 
defect type (and m) 


Status 


XY/RP 1 /0(2) 
& step model 


2 


Z, vortices (m = 0) 


BKT transition (see section |H| 


0(n) 
for n > 3 


2 


No defects 


No transition (majority opinion - e.g. 

[muni mug, but see also [in no ma ezd 


RP™- 1 
for n > 3 


2 


Z2, vortices V n 
(m = 0) 


No transition [|S] (perturbation theory); 
l st -order transition [ I^Hl IH03 ( n ~ ¥ °°); 
BKT or 2 nd -order transition 

for n = 3 [EZHS3EHE21; 

No transition for n = 4 [I65j 


U[n) 


Q 
O 


/t — l. lu^ voi Lices 
(m = 1); n = 3: Z, 
monopoles (m = 0); 
no defects for n > 4 


z -oiuei LianbiLKjiis 

(see [H2] and references therein) 


RP™- 1 
for n > 3 


3 


n = 3 only: Z, 
monopoles (m = 0); 
n > 3: Z2, vortices/ 
disclinations (m = 1) 


2 nd -order transition [[^ (from 
perturbation theory); 
l st -order transition [56 (n — > 00); 
l st -order transition [EHl E21 ESI (" = 3); 
1 -order transition [1541 155j (n = 4) 



was demonstrated that the apparent phase transition may be completely eliminated. 
The suppression is achieved by the introduction of a chemical potential term asso- 
ciated with the defects, making the formation of topological charges energetically 
expensive. 

The suppression of the monopoles in the three-dimensional Heisenberg model (which 
is also non-Abelian) also leads to the disappearance of the transition there [EH]- Sim- 
ilar work was applied to the nematic/isotropic transition of the three-dimensional 
RP 2 model in [ E2] • The first-order transition between the nematic and isotropic 
phases weakens as the disclination core energy is increased and eventually the tran- 
sition line splits into two continuous ones separating three distinct phases with a 
new topologically-ordered isotropic phase between the nematic and conventionally 
isotropic ones. 
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However, while the work of [ E3 EH1 EH E2j favours a second-order or BKT-type 
transition with divergent correlation length in the n = 3 case, and (see also [ I64j ) . 
another study [ 153] of the RP n_1 defects in a d = 2 model of tops favours the absence 
of a true phase transition there (at least for n — 4). Instead it is claimed that there 
is a crossover in the correlation length. There it is argued that the defects disorder 
the system for all temperatures and the correlation length remains finite. 

Therefore, the situation in the d = 2 RP n_1 models is still not satisfactorily 
clear and the precise nature of the phase transitions in these models is still under 
question. 

In table EJ a summary of the status of the standard 0{n) and RP n_1 models and 
the associated homotopy groups is given, together with a selection of recent papers. 

9. Highly nonlinear models 

Universality is the notion that the existence and type of phase transition in a 
model, and the critical exponents that describe it, depend only on the dimension, 
the symmetries present and the range of interaction. One can broaden the scope the 
models discussed herein, while maintaining these three characteristics, by altering 
the interaction to the form 

E = -Y,Pk(s l s J ), (18) 

where Pk{x) is the k th Legendre polynomial. Then k = 1 and k = 2 correspond to the 
0(n) and RP n ~ l models appropriately. The k = 4 version has P±{x) cx 35x 4 — 30s 2 +3 
and has the same symmetries as the k = 2 version but a higher degree of nonlinearity. 

In d = 3 dimensions, the inclusion of the P4 interaction enhances the first-order 
transition present in the k = 2 model for n = 3 (i.e., the RP 2 model) [ 66:. Similarly, 
in two dimensions, the n = 3, P4 model was found to have a strong first-order 
transition [EZj- In two dimensions, the claimed (but disputed) continuous transition 
in the n = 3, P2 model (i.e., the RP 2 model) and the first-order transition in the 
n = 3, P4 model can be eliminated by suppression of defects [152]. Finally, in two 
dimensions, a BKT-like transition was recently observed in the n = 2, P4 model in 
[ 68J. There, the crucial role of topology in such models is again exhibited and new 
questions are raised regarding the detailed nature of that transition. These latest 
results demonstrate yet again the importance and ubiquity of topologically mediated 
transitions in contemporary statistical mechanics. 

Recently another type of non-standard n-vector model has come into vogue. This 
is given by a nonlinear Hamiltonian of the type 

P = -£(1 + ^) P . (19) 

M 

If p = 1, this recovers the standard 0(n) models. There, for d = 2, the n = 2 
version has the BKT transition discussed above, while for n > 2 there is no positive- 
temperature transition. In [HSU, [EH and [E] numerical evidence for the existence 
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Table 6. Summary of the status of some highly (i.e., sufficiently) nonlinear models 
and some recent developments considered in this paper. Here m is the dimension 
of the defects. 



Model 


d 


Homotopy group, 
defect type (and m) 


Status 


n = 2, P 4 


2 


Z, vortices (m = 0) 


BKT-like transition [ 168] 


n = 3, Pi 


2 


Z2, vortices (m = 0) 


1 st . order transition [1621 167j 


n = 3, P4 


3 


Z, monopoles (m = 0) 
Z2, vortices (m = 1) 


1 -order transition [I66j 


Nonlinear 0(n) 


2 


n = 2: Z, vortices (m = 0); 
no defects if n > 3 


l st -order transition for 

n = 2,3,00 [Hanoi HQ 

and generally [ E3 ESI 


Nonlinear RP"" 1 


2 


^2) vortices (m = 0) 


1 -order transition [ I73j 


Nonlinear 0{n) 


3 


n = 2: Z, vortices (m = 1); 
n = 3: Z, monopoles (m = 0); 
no defects for n > 4 


1 -order transition [I72j 


Nonlinear KP n ~ l 


3 


n = 3 only: Z, monopoles (m = 0); 
n > 3: Z2, vortices (m = 1) 


l st -order transition [I73j 



of first-order transitions in the n = 2, n = 3 and n = 00 (spherical) models respec- 
tively in d = 2 dimensions for sufficiently large p was proffered. A rigorous proof 
of the existence of first-order transitions in these models was given by van Enter 
and Shlosman [ E2] These transitions can occur in d = 2 dimensions despite the 
implications for zero magnetization coming from the Mermin- Wagner theorem [ |3] 
and despite the high-p models sharing the same dimension, symmetries range of 
interaction as the standard p = 1 versions. 

In three dimensions, where the Mermin- Wagner theorem does not apply, the 
standard 0(ri) models, given by (fTU|) with p — 1, have second-order transitions 
described by n-dependent critical exponents. Again, the traditional notions of uni- 
versality would imply that such behaviour is independent of p. However, sufficient 
nonlinearity (sufficiently large p) can cause first-order transitions in these models 
and the rigorous proof of [ E2] extends to these d = 3 cases too. 

In [ES], it i s rigorously shown that various sufficiently nonlinear models of the 
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RP n 1 type also exhibit first-order transitions in d > 2 dimensions. Here, the Hamil- 
tonian is of the form 



From these recent developments, it is clear that the notion of universality has to be 
extended to include the degree of nonlinearity as a significant factor. For an overview 
of recent developments in nonlinear models, see [[71] and table El 

10. Conclusions 

Topology plays an important role in condensed matter physics. In particular, 
homotopy theory facilitates the understanding and classification of the conditions 
which permit the existence of topological defects - domain walls, vortices, monopoles 
and so on. After an expeditious introduction to the essentials of homotopy theory, a 
variety of topologically mediated phase transitions has been surveyed, starting with 
the BKT transition in the XY model in two dimensions - the paradigm for such 
studies. After three decades of work, the detailed perturbative RG predictions of 
[ 01 E] have been confirmed in that model through taking multiplicative logarithic 
corrections into account. Furthermore, recent progress in confirming the vortex- 
binding mechanism as that mediating the phase transition in the XY model has 
been reviewed to complete an account of the current status of one of the most 
beautiful and remarkable models of theoretical physics. 

Besides the XY model, a number of other phase transitions have been briefly 
examined, some of which are topologically-mediated. Some open problems have been 
highlighted, the resolution of which now stands at the forefront of modern analytical 
and numerical investigations into statistical mechanics. 
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